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Data Assimilation and Transport Modeling in Terrestrial and Planetary 

Atmospheres 


Data assimilation is a blanket term used to describe a number of techniques for retrieving 
important physical parameters from observational data, subject to constraints imposed by prior 
knowledge (such as, in the case of meteorology, the primitive equations that govern 
atmospheric motion). Since these newly developed methods make efficient use of 
computational resources, they are of great importance in the interpretation of the voluminous 
datasets that are now produced by satellite missions. 

As proposed, these techniques have been applied to the study of the Martian and terrestrial 
atmospheres based on available satellite observations. In addition, a sophisticated 
hydrodynamic model (non-hydrostatic, and therefore applicable to the study of the interiors of 
the giant planets) has also been developed and successfully applied to the study of tidally 
induced motions in Jupiter. 

Mars 

For study of the huge meteorological dataset returned by the Mars Global Surveyor (MGS), a 
new baroclinic spectral general circulation model (GCM) has been utilized. It includes well- 
understood dynamics and the ability to be truncated to the resolution of the MGS observations 
|4). Given the available observations for the mapping phase of the MGS mission, the model 
state is thus overdetermined by the observations [2|. Good agreement between the analyzed 
state and the Thermal Emission Spectrometer (TES) retrievals is obtained (rms error ~ 2K; 0 
bias) 1 1 1 1. The straightforward analysis also provides a satisfactory forecast (rms error < 4K) 
112), but a future sequential assimilation will provide even better forecasts [14] and additional 
information about the sensitivities of the Martian atmospheric state (that should lead to 
improved future missions) [17J. It is now possible to assimilate the observed TES radiances 
directly 1 19), though some systematic differences between the individual TES retrievals and the 
assimilation analysis are evident [21 ]. A full year's analysis of TES-derived meteorology 
shows generally good agreement with traditional GCM calculations [22]: the average 
temperatures are generally warmer by about 10 K, the equator-pole temperature gradient and 
thus the mean zonal wind are consistent with model predictions (and with the gradient wind 
approximation). The meridional Hadley cell is about as predicted, but the transition from 
northern summer to southern summer circulation is extremely rapid. With the goal of 
validating this TES dataset against other meteorological measurements, a systematic 
comparison with radio occultation profiles has been conducted [24|. Relatively large rms 
errors (~ 6K) are an indication that the near-surface TES retrievals and the assimilation model 
ground temperature and boundary layer forward models should be revised considerably. In 
particular, low-level inversions that are common in the (nighttime) occultation profiles are not 
found in the analyses. 
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The assimilation of TES-derived water vapor column amounts awaits improved Martian 
atmospheric transport models. Current models are able to produce a closed water cycle that 
resembles Mars Atmospheric Water Detector (MAWD) and TES observations [5], However, 
these models show very small diurnal variations of water vapor column [1 J and this lack of 
variation cannot be verified from the limited TES observing pattern. Nevertheless, the model 
can be used to simulate near-surface ground water (which exists as adsorbed water, vapor, and 
ice). The computed amounts are smaller than what has generally been derived from models 
which ignore the seasonal variation of atmospheric column [3]. The Mars Odyssey (MO) 
promises to provide the observations that can verify the new predictions [16J of smaller near- 
surface water content in the southern hemisphere than in the northern hemisphere. (As of 
spring 2002, however, MO has not yet seen below the seasonal frost deposits in the northern 
hemisphere.) 

Earth 

Data assimilation of terrestrial aerosol observations from space is difficult because of the many 
variables associated with the size distribution. In an effort to simplify this process, a moment- 
oriented approach to the modeling of stratospheric aerosol sedimentation was developed. This 
showed that the timescale for evolution of volcanically injected sulfuric acid aerosol clouds in 
the stratosphere lengthens with time [6], consistent with observations of visible extinction by 
SAGE II and at infrared wavelengths by HALOE [7]. The predictions of power-law decay of 
extinction are borne out by the data [8|. This in turn allows the formulation of scaling laws 
that predict the integrated climatic impact of a volcanic eruption as a function of total eruptive 
mass |9|. A detailed study of the ratios of infrared to visible extinction, however, shows much 
less time variation than predicted 1 10), consistent with an effective radius of the post-Pinatubo 
volcanic aerosol much smaller than expected [15J. This has spurred the development of a full 
assimilation model for the study of satellite measurements of stratospheric aerosol (20). The 
model analyzes the full range of the particle size distribution, the mean meridional Brewer- 
Dobson circulation in the stratosphere, and source and sink terms of sulfur dioxide on a 
monthly mean basis [23]. 

Jupiter 

The innermost Galilean satellites of Jupiter all show interesting phenomena that can ultimately 
be traced to anomalously high heating rates: the volcanoes of Io; a probable liquid water layer 
in Europa; and an active magnetic dynamo in Ganymede. The ultimate source of this energy is 
tidal dissipation in Jupiter: Io raises a bulge on the fluid planet (analogous to the ocean tides 
raised by our moon) which exerts a torque back on Io, accelerating it in its orbital motion, 
leading to strong interactions with the other satellites. This torque is only possible if some 
(friction-like) dissipation mechanism delays the occurrence of the bulge until after the passage 
of the exciting satellite. No adequate source of this dissipation in Jupiter has yet been 
identified, largely because of the lack of a complete hydrodynamic model of this planet [13). 
Such a model has now been developed and a number of previously troubling questions about 
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Jovian tidal dynamics have now been answered f 18]. Among these are the correct boundary 
condition that connects the adiabatic envelope of the planet to its stably stratified atmosphere; 
the role of a stable layer in the envelope in amplifying tidal motions; the correct definition of 
the tidal dissipation factor (Q) in a compressible body; and the relation of the atmospheric 
wave flux to the tidal torque. Still under investigation is the question of where the angular 
momentum exchanged in the tidal interaction is deposited in the planet. 
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Appendix 


Subject Inventions Certification 

There were no subject inventions required to be disclosed to NASA which resulted from 
this work. There were no subcontracts awarded under this Cooperative Agreement. 



